Maintaining Quantum Coherence in the Presence of Noise by State Monitoring 
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Unsharp POVM measurements allow the estimation and tracking of quantum wavefunctions in 
real-time with minimal disruption of the dynamics. Here we demonstrate that high fidelity state 
monitoring, and hence quantum control, is possible even in the presence of classical dephasing 
and amplitude noise, by simulating such measurements on a two-level system undergoing Rabi 
oscillations. Finite estimation fidelity is found to persist indefinitely long after the decoherence 
times set by the noise fields in the absence of measurement. 



o 

(N 

> 
O 



o 

oo 



X 



Maintaining high-fidelity quantum control is a central 
requirement in a variety of technologies ranging from 
nuclear magnetic resonance to quantum based precision 
measurement 0]. Quantum control is usually restricted 
to a finite time window as a result of the unavoidable in- 
fluence of decohering environments, and the control life- 
time is often extended through the use of decoherence 
free subspaces 0] , or dynamical decoupling 0] . 

In this communication we discuss a scheme which relies 
on a sequence of consecutive POVM (Positive Operator- 
Valued Measure) measurements to maintain quantum 
control. We demonstrate that the time-evolution of a 
driven, isolated two-level quantum system, subject to 
classical dephasing and amplitude noise, can be moni- 
tored long beyond its Rabi coherence time. In fact, the 
wavefunction can in principle be tracked indefinitely with 
finite fidelity, unlike systems controlled by dynamical de- 
coupling that ultimately undergo complete loss of coher- 
ence. The control scheme relies on periodic application 
of special POVM measurements said to be "unsharp" (3) ■ 
Such measurements have previously been shown to allow 
faithful monitoring of Rabi oscillations if the general form 
of a time- independent Hamiltonian is known [^, Q and 
no external noise is present. 

A scheme for updating a state estimate during contin- 
uous measurements Q - the continuum limit of the tech- 
nique employed here - has been presented in In that 
scheme the state evolution of the system, its estimated 
state, and the measurement readout is described by three 
coupled stochastic differential equations, which indicate 
that the state estimate converges to the real state for a 
broad class of systems, cf. @. Continuous measurements 
have also been shown to drive statistical mixtures of spa- 
tial wavepackets into pure states, which can be entirely 
determined by the measurement record alone (Io| . 

Specific experimental implementations of unsharp 
measurements have been suggested in the context of 



Bose-Einstein Condensates cavity QED [12| and 
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coupled quantum dots 13, [Hj]. Several realizations of 
the related topic of "weak- value" measurement have been 
demonstrated through measurements of photon momen- 
tum [Il-Gll. In addition, experiments of "continuous 
weak measurement" were implemented using a cold ce- 
sium vapor [nH - These realizations all employed ensem- 
ble measurements, while we here show monitoring of a 
single, isolated quantum system by repeated measure- 
ment, as the system evolves. 

We consider a two-level system undergoing Rabi- 
oscillations. In a frame rotating at the two-level tran- 
sition frequency it evolves under the Hamiltonian 



Hr 



H — CTa. 



(1) 



where a x is the Pauli matrix which generates rotations 
about the x-axis and £Ir the Rabi frequency, which is 
assumed to be known. At the same time we assume the 
system is under the influence of random classical noise 
fields p(t) and a(t), causing dephasing and amplitude 
fluctuations, respectively, through a noise Hamiltonian 



H N = h/3(t)a z + Ha(t)a x 



(2) 



Each noise field is characterized by a power spectrum 
which is related to its autocorrelation function 



through 



P({u) = J C (2) (r)e l " 



r)dt 



r dr 



(3) 



(4) 



where = a(t),0(t). 

The estimation strategy rests on carrying out POVM 
measurements periodically and updating the state 
estimate based on the measurement outcomes. Quite 
generally, a POVM measurement with outcome n, which 
was carried out on a system in the state \ip), will result 
in a state after the measurement given by 



(5) 
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Here M n is the so-called "Kraus operator" corresponding 
to the measurement outcome n, and 



p(n|V) = (^MlMnW 



(6) 



is the probability to detect outcome n, conditioned on 
the system being in state 

In an estimation experiment a sequence of periodic 
measurements, with period r, are applied to the system 
as it evolves in time [6|. Despite the dynamics, the state 
change due to the measurement can still be described by 
Eq. if each measurement is executed much faster than 
all other dynamical timescales (impulsive measurement 
approximation). In between measurements the time evo- 
lution is described by the operator 



Uh = T 



exp(- 



tj+T 



(H R + H N (t))dt) 



(7) 



where T is the time-ordering operator. At t = Nt, after 
N measurements, the system is, up to the appropriate 
normalization constant, in the state 

IV(JVr)) = M nN U N M nN _ 1 U N - 1 ...M ni U 1 \^)- (8) 

To estimate the state of the system the same sequence 
of operators corresponding to the measured outcomes in 
Eq. (|8) are applied to an initial guess \ip e ) (cp. 0), but 

(1) the initial estimate of the state \ip e ) can be taken as 
an arbitrary state vector on the Bloch sphere and 

(2) in between measurements the state estimate is as- 
sumed to evolve only through the Hamiltonian Eq. (flj, 
since the experimenter does not know what the instanta- 
neous values of the noise fields are. As we'll see in what 
follows it is still possible to estimate the state of the sys- 
tem without detailed knowledge of the noise fields. Our 
approach differs from [H[ where the Rabi frequency was 
assumed to be unknown and one of the aims was to de- 
termine its value through a Bayesian estimator in the 
absence of noise. 

We now define two projectors P + = i(l + f • a) 
and P- = i(l — f • <r), where 1 is the identity oper- 
ator, f = (6, C, x) a unit vector on the Bloch sphere, 
and a = (<r x , a y , a z ). In terms of these it is possi- 
ble to construct POVM measurement operators Mq = 
y/pEP+ + VT - J* and Ah = VT~^ P+ + ^/pE P-, 
related via M^M + m\m x = 1, and < p < 0.5. 
The strength of a single measurement is quantified by 
Ap = (1 — po) —po = l — 2po [2l|. However, the strength 
of a sequence of measurements depends also on the period 
t between two consecutive measurements. For fixed Ap 
a shorter (longer) period r means a stronger (weaker) 
influence of the sequential measurement. The strength 
of the state disturbance due to this sequential measure- 
ment is best quantified by the rate j m = 1/ r TO with 
r m = 2r/(Ap) 2 [2l|. The strength 7 m is the expected 
rate at which an arbitrary initial state is reduced to an 
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FIG. 1: Wavefunction estimation in the absence of noise for a 
specific realization of a(t) and /3(t) (single run), (a) Expecta- 
tion value (a z ) for blackline - true expectation value, red line 
- estimated expectation value, dashed blue line - expectation 
value in the absence of measurements, (b) Estimation fidelity. 



eigenstate of the measured observable, in the absence of 
dynamics other than measurement 0. 

To set the stage we illustrate the method in the absence 
of noise, i.e. /3(t) = a(t) = 0. We simulate an experi- 
ment in which we choose Ap = 0.2 and f = (0, 0, 1) which 
corresponds to an unsharp measurement of er 2 (cp. [H[). 
We carry out a measurement every r = Tr/10, where 
Tr = is the Rabi-period. The resulting measure- 

ment strength is thus smaller than the Rabi-frequency 
(7 m = f2fl/107r), which is required in order not to dis- 
turb the oscillations too strongly. For a measurement 
strength 7 m much greater than the state would be 
projected onto an eigenstate of the observable a z before a 
single oscillation can take place, and the dynamics would 
freeze (similarly to the Quantum Zcno Effect 



22(1). The 

result of each measurement is chosen at random, com- 
mensurately with the probabilities prescribed by Eq. ^ . 
The initial state estimate is chosen orthogonal to the ini- 
tial state vector, a limiting case for which the estimation 
procedure might be expected to have some difficulty 

In Fig. QJa) we plot the expectation value of (a z ) for 
the true state (black line) and the state estimate (red 
line) for one single run of the measurement experiment 
(i.e. one specific realization of a(t) and f3(t)). The state 
(black line) undergoes Rabi oscillations, but with mea- 
surement induced random phase shifts as compared to 
the undisturbed oscillations (dashed blue line). The os- 
cillations including the influence of the measurement are 
monitored accurately by the estimate (red line) after 
about 6 Rabi periods. After this time not only the ex- 
pectation value of the measured observable a z with re- 
spect to the true and estimated state, but also the states 
themselves coincide as is indicated by the plot of the es- 
timation fidelity F(t) = \(ip e st\ip}\ 2 in Fig.QJb). Asymp- 
totically the fidelity tends to unity, indicating the perfect 
state monitoring of a single system, in real time, in the 
absence of noise. 
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FIG. 2: State estimation in the presence of dephasing and am- 
plitude noise, (a) Expectation value of a x and (b) expectation 
value of <t z for a specific realization of a(t) and P(t). Here we 
used A/3 = 0.05, Aa = 0.005, Ap = 0.2, and f = (0.43, 0, 0.9). 
Note (a) and (b) have different time axes. 
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FIG. 3: State estimation in the presence of classical noise, (a) 
Red line - Estimation fidelity for a single run (corresponding 
to Fig. [2| , Black line - expected fidelity obtained by averaging 
f 000 runs, (b) Blue line - Rabi oscillations showing monotonic 
loss of coherence in the absence of unsharp measurements. 



Now consider a more realistic situation in which the 
two- level system is not isolated, but subject to random 
classical noise as described by the Hamiltonian Eq. ([2]). 
As an example we assume the noise fields both have a 
power spectrum P^(cj) = A^/lj, since this "one-over-f 
noise is ubiquitous in many systems. For concreteness we 
choose a lower cutoff of w = 0.01 and a high-frequency 
cutoff of uj = 10, where the frequency is specified in units 
of fi/j. In accordance with Eqs. §5$ and we generate a 
specific noise trajectory by summing over different spec- 
tral components, weighing each with the square root of 
the noise power: £(t) = v^7(^i) cos + fa). Each 

i 

spectral component contains a random phase factor fa, 
allowed to vary between [0, 2tt], and assumed to be delta 
correlated. Each noise trajectory, a(t) and /3(t), is nor- 
malized so that their root-mean-square deviations are re- 



spectively one hundredth and one tenth of the drive field 
amplitude, Aa = 0.005 and A/3 = 0.05. We use the 
same measurement strength as before, but an observable 
f • a with a finite x- component: f = (0.43,0,0.9), since 
the noise is expected to tip the Bloch vector out of the 
yz-plane. 

Figure [2] displays the evolution of the expectation val- 
ues (a) (a x ) and (b) (a z ) for the true state (black lines) 
and the estimate (red lines), again for a single run of 
the experiment. The amplitude of the Rabi oscillations, 
Fig-HJb), although modulated by the noise, does not de- 
crease permanently and the estimate succeeds in tracking 
both components. The red line in Fig.[3][a) shows the es- 
timation fidelity corresponding to the single run of the 
experiment of Fig. ([2]). Despite the noise the state es- 
timate quickly approaches the real state, although the 
fidelity does not converge completely to unity. Instead, 
it exhibits random excursions away from unity, which at 
long times are centered around an average, asymptotic 
value. To find this value we execute 1000 runs of the 
experiment with the same initial conditions and average 
over the resulting fidelities, leading to the black curve 
in Fig. [UJa). We've found empirically that this average 
fidelity F is well described by 

F(t) =F (l-e-*/ TE ). (9) 

In Eq. (|9|), Fq is the asymptotic estimation fidelity and te 
the estimation time. By fitting Eq. ^ to the simulated 
result we extract an estimation time of te = 3.7Tr and 
an asymptotic fidelity of F = 0.98. For comparison the 
blue curve in Fig. E^b) plots the result of an average over 
1000 simulated runs of the experiment in the absence 
of measurements, but with the same noise source as in 
(a). It shows the decay of Rabi oscillations to about 
half full amplitude over the same time span due to the 
noise. Labelling tr the characteristic decay time of Rabi 
oscillations, we remark that Eq. ((9]) holds accurately only 
when r m <C tr. For r m > tr the asymptotic approach is 
no longer simply exponential. 

The results of Figs. (a) and (b) taken together im- 
ply that in any single run of the experiment the state 
can be estimated at all times after convergence with an 
average of 98% fidelity, by a pure state \if> e (t)). On the 
other hand, in the absence of measurements the state 
would lose coherence due to the noise, and evolve into 
a statistical mixture as evidenced by the decay of Rabi 
oscillations of the ensemble average shown in Fig. |3] (b). 
This constitutes the main result of this letter. The fi- 
delity can be operationally tested at the end of a run, by 
deducing from the state estimate the appropriate unitary 
rotations needed to place the system in the state | f) , say, 
where it will then be detected with 98% probability. As 
such, the experimenter has maintained quantum control 
by monitoring of the state evolution, despite the noise. 

Finally we study the effectiveness of the estimation 
process as a function of the measurement strength, j m , 
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FIG. 4: Effect of measurement strength on state estimation in environments with different noise strengths: diamonds - No 
noise, circies - A/3 = 0.05, Aa = 0.005, squares A/3 = 0.1, Aa = 0.01. (a) Asymptotic fidelity, Fo, as a function of measurement 
strength, (b) Convergence time, te, as a function of measurement strength. For the three weakest 7 m 's on the curves with 
noise, we plot in (a) the asymptotic fidelity and in (b) that time which it takes to reach a fraction 1 — e -1 of that fidelity, even 
though the corresponding fidelity curves aren't strictly exponential. 



when noise is present. The simulations are repeated for 
different values of j m , and still choosing f = (0.43, 0, 0.9) 
in each case. Figures Ufa) and (b) respectively plot the 
estimation fidelities and convergence times as a func- 
tion of 7 m for different noise strengths: diamonds - 
no noise, circles - A/3 = 0.05, Aa = 0.005, squares - 
A/3 = 0.1, Aa = 0.01. When noise is present, the asymp- 
totic fidelity monotonically decreases as the measurement 
strength becomes weaker, and approaches Fo = 0.5 as 
7 m — > 0. This is consistent with the average fidelity 
obtained when taking random guesses for the state es- 
timate. Simultaneously, the convergence time increases 
as the measurement becomes weaker, but plateaus to a 
finite value as 7 m — > 0. By contrast, in the absence of 
noise the fidelity always approaches Fq = 1, but the con- 
vergence time increases indefinitely as the measurement 
strength weakens, as can be seen in Fig.|4jb). 

The trends observed in Fig. [4] emphasize that the ap- 
propriate timescales need to be obeyed for the measure- 
ment scheme to work. The sequential measurement must 
be weak enough not to freeze the dynamics, but strong 
enough to enable a high fidelity estimate before the noise 
randomizes the system, i.e. Tr <C r m <C tr. 

In conclusion, we remark that this study was carried 
out in a regime of comparatively strong noise, namely 
Hr/A/3 = 5 ~ 10. With stronger drive fields higher 
asymptotic fidelities can be expected for the same mea- 
surement strengths considered here. For example, we 
find that with CLr/A/3 = 100, ft R /Aa = 1000, Ap = 0.1, 
r = (0.43,0,0.9) ,that F = 0.999 and r = 15.8Tr. It is 
encouraging that the estimation procedure described here 
predicts finite estimation fidelity despite the presence of 
random classical noise. This opens the way for quantum 
control techniques that monitor wavefunction dynamics 
beyond the limitations set by decoherence processes in 
the absence of unsharp measurements. 
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